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Abstract 

We consider a disordered asymmetric exclusion process in which randomly chosen sites do not 
conserve particle number. The model is motivated by features of many interacting molecular 
motors such as RNA polymerases. We solve the steady state exactly in the two limits of infinite 
and vanishing non-conserving rates. The first limit is used as an approximation to large but finite 
rates and allows the study of Griffiths singularities in a nonequilibrium steady state despite the 
absence of any transition in the pure model. The disorder is also shown to induce a stretched 
exponential decay of system density with stretching exponent = 2/5. 
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I. INTRODUCTION 



Driven diffusive systems serve as simple models for collective phenomena ranging from 
traffic flow to molecular motors. Moreover, they provide tractable examples of systems far 
from thermal equilibrium. Studies of one- dimensional driven diffusive systems have shown 
that many interesting phenomena, which are typically not observed in one dimensional 
systems in thermal equilibrium, exist. Prominent examples are boundary induced phase 
transitions and spontaneous symmetry breaking, for reviews see e.g. 3, [J. 

Most studies have considered systems in which the dynamics are the same everywhere in 
the system or systems where the dynamics are modified only at the boundaries. However, 
when trying to relate these systems to many interacting molecular motors, the effects of non- 
conservation and disorder (i.e. spatial heterogeneity in the dynamics) can not be ignored in 
many cases. 

Indeed there have been some studies on the effects of disorder on driven diffusive systems. 
For example, the effect of assigning a disordered quenched rate to each particle was studied 
in J, 0, Ifi , 21 Isl El on a ring geometry. Exact solutions show that at high enough densities 
a macroscopic number of particles jam behind the slowest particle in the system. The 
phase transition between the jammed and non-jammed phase is similar to a Bose-Einstein 
condensation. Work has also been done on an asymmetric exclusion process on a_ ring 
where the quenched hopping rates between neighboring sites are drawn at random 
For molecu 
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ar motors moving along a disordered substrate this seems to be the relevant 
ll^ . It was argued, based on numerics and mean- field solutions, that at high 



densities the system phase separates into a region of high density coexisting with a low 
density region. Finally, the combined effect of random hopping rates and open boundary 
conditions was considered in it was argued using numerics that the location of 



phase transition lines may be sample dependent. In [1^ mean-field arguments and numerics 
indicate the existence of shifts in phase boundaries which, by analogy with equilibrium 
systems, are expected to be accompanied by emergent Griffiths regions. A review of the 
effects of disorder in exclusion models has been given in jis| . 

In this paper we consider another type of disorder. We study an asymmetric exclusion 
process (ASEP) where non- conserving sites are chosen at random along the lattice. At 
these sites particles may attach and detach with specified rates which may also be drawn 
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at random. Thus there are two components to the disorder. A feature of this disorder 
is that it allows a detailed account of the way in which Griffiths singularities can arise in 
nonequilibrium steady states. In equilibrium the mechanism leading to Griffiths singularities 
is well understood: the disorder, e.g. dilution, breaks the system into pure regions and large 
pure regions may give rise to the exponentially suppressed Griffiths singularities. In the 
present the case, in the limit of high attachment and detachment rates, the non- conserving 
sites break the system into driven conserving domains. 

Non-conservation of particles in driven systems without disorder has previously been 
considered in the context of molecular motors. The idea is that molecular motors move in a 
preferred direction along a filament and are able to attach and detach from the filament. In 



the works so far all sites are non-conserving 



16, 
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1^^^. The motivation for the model 



we study here comes from the fact that some molecular motors only attach and detach at 
certain sites. 

More specifically, we give a very simplified description of many interacting RNA poly- 
merase (RNAp) motors acting on a prokaryotic DNA in vitro. Prokaryotic RNA polymerase 
can initiate without regulatory proteins. Namely, RNAp left in a solution with DNA can 
produce RNA even if the specific protein which regulates its action (for example, by enhanc- 
ing or reducing the initiation rates) is present. They can enter and leave the DNA in order 
to transcribe RNA molecules at specific sites, referred to as promotor and termination sites 
respectively. In the language of the lattice model we consider a binding of a RNAp to a pro- 
motor corresponds to a particle entering the system. The unbinding at the termination site 
corresponds to a particle leaving the system. In the absence of regulatory proteins the rate 
of entering the DNA depends on the details of the promotor sites. In such systems the RNA 
polymerase do not usually move from one gene to another. In the lattice model this would 
correspond to particles not moving from one stretch of conserving sites to a neighboring 
one i.e. the limit in which the detachment rates are large at the non-conserving sites. We 
comment that in principle RNA polymerase may move in different directions along the DNA 
when transcribing different messenger RNAs, corresponding to particles moving in different 
directions along different conserving stretches. However, in the limit when the detachment 
rate is large this will not influence most of the results described in the paper. Of course, 
our assumption of randomly distributed lengths of genes (or conserving segments) is not 
expected to hold. However, the model provides a starting point for analyzing more realistic 
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situations. 

The paper is organised as follows. In Section |Hj we define the model and discuss two lim- 
its which are exactly soluble. In Section UTTl we show that the disorder induces Griffiths like 
singularities as the rates for entering and leaving the lattice are changed. More significantly, 
it is shown in Section HV1 that the presence of the non-conserving sites leads to anomalous 
relaxation of the system toward the steady state. Specifically, we argue that decay of mea- 
surable quantities decay as a function of time, t, as a stretched exponential exp(— ct^), where 
c is a non-universal constant and = 2/5. The results are verified numerically. We conclude 
in Section fVl 

II. MODEL 

The model we consider is a disordered generalisation of the ASEP. The pure ASEP is 
defined on a one-dimensional lattice containing L sites and with periodic boundaries. The 
lattice is occupied by particles subject to an exclusion interaction, which prohibits multiple 
occupancy of any site. These particles hop with rate one to the nearest neighbour site to the 
right, provided it is empty, and so the total particle number N is conserved. We introduce 
non-conservation into this model by allowing, at certain sites (which we will call 'disorder' 
sites), processes which do not conserve the total particle number N. Hence each site I 
(I = 1, . . . , L) in the pure model remains a pure site with probability p, or becomes a disorder 
site with probability (1 — p). Now, at the disorder sites, labelled by j = 1, . . . , P, particles 
attach with rate Cj or detach with rate aj. In general, we wish to consider heterogeneous 
rates for the non- conserving processes. 

To study the model we first consider limits which can be solved exactly. Later, using 
numerics, we argue that the results are generic. Exact solubility arises when the steady 
state densities at the disorder sites are determined solely by the attachment and detachment 
processes. In these cases, since the system is composed of conserving domains of the chain in 
contact with disorder sites at the boundaries of each domain, the steady state can be written 
as a product of boundary driven ASEPs in which the densities of the boundary reservoirs 
are given by the disorder site densities pj. 

Before turning to the disordered ASEP under consideration we recap some facts, which 
will be useful later, about the boundary driven ASEP. For the boundary driven ASEP, in 
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which particles are injected at the left-hand boundary site with rate a (provided it is vacant) 
and removed from the right-hand boundary site with rate (3, exact steady state weights for 
particle configurations can be obtained using a matrix product ansatz 20]. In this ansatz, 
particle configurations are represented as a product of matrices Xi ■ ■ ■ where Xi = D 
(E) if site I is occupied (vacant). The steady state weight of a configuration is given by 
(a\X\ - ■ ■ Xl\(5) provided the matrices D and E and the vectors (a| and \(3) satisfy the 
relations 

DE = D + E = C, a(a\E=(a\ and /3D\/3) = \/3) . (1) 

From these relations exact expressions for the normalisation (a\C L \/3) can be derived which 
show that the model undergoes a second order phase transition: when both a and (3 > 1/2 
the system is in a maximum current phase, otherwise it is in one of two low current phases. 
The phase transition between the low current phases is first order. The rates a and 1 — (3 
represent the densities of particles in reservoirs connected to the boundary sites. 

Next, we use the known results for the boundary driven ASEP to study the disordered 
case. As stated above, there are limits where the steady state weight of the disordered model 
factorizes into a product over boundary driven ASEP weights. The two exactly soluble limits 
are: 

• Cj, aj — > oo, with Cj/cij fixed. 

In this limit, each disorder site j acquires a density pj determined solely by Cj and aj 
which obeys the equation of motion, 

dpj 

= Cj {l - Pj ) + a jPj . (2) 

Therefore in the steady state, 

Cj + aj 

If we define rij to be the number of sites between disorder sites j and j + 1 (i.e. 
the length of the j-th conserving domain), then the normalisation ZL({rij}) (which is 
the sum over the steady state weights of all particle configurations on sites excluding 
the disorder sites), for a given configuration of the disorder sites {fij} = rii, . . . ,np, 
factorises into a product over normalisations for the boundary driven ASEP: 

p 

Z L ({n J }) = l[(p,\C^\l- Pj+1 ) (4) 



and where pp+i = p\. This is the relevant limit for the model molecular motors 
discussed in the introduction. 



Cj, Clj 



0, with Cj/cij fixed. 

In this limit the time between each attachment/detachment event tends to infinity. 
Therefore, after each event the system reaches a homogeneous steady state of the 
pure ASEP with periodic boundaries. Thus the system density p = N/L, satisfies the 
equation of motion 

dp 



dt 



£ 

.3=1 



(1-P) + 



j=i 



p 



Therefore in the steady state the system density is given by 



P 



E 



j=i c j 



E 7 =l( C j + a 3 



(5) 



(6) 



Because the steady state is homogeneous, all sites, including disorder sites, have the 
same steady state density p. Moreover, the steady factorises and there are no correla- 
tions between sites. Therefore, one can still write the normalisation in a form similar 
to ®: 

p 

^(K})=n<H^ii-p>. (7) 

3=1 

This is because in this case D and E are given by the scalars 1/(1 — p) and 1/p. 

In the following we will consider the first limit, Cj, aj — > oo. We use the factorised form 
with pj given by Q as an approximation for the case where large but finite 

which is relevant for the model of molecular motors. This approximation has a mean-field 
character, in the sense that correlations are factorised about the disorder sites, however all 
correlations within conserving domains are retained. 

III. GRIFFITHS SINGULARITIES 



We can exploit known properties of the normalisation of the boundary driven ASEP to 



demonstrate the existence of Griffiths-type singularities in the disordered ASEP 
an illustrative example, we consider binary disorder at disorder sites, such that 

u with probability q, 
v with probability 1 — q. 
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221. As 



P.i 



aj + Cj 



This is the simplest choice of disorder for which Griffiths singularities occur. Generalizations 
to more complicated situations are straightforward. 
Using (JD) the steady state normalization satisfies 

p 

]nZ L = ^2W nj { Pji l- Pj+1 ) , (8) 

i=i 

where W n (pj, 1 — Pj+i) = ln(/?j|C nj '|l — Pj+i). In order to perform the disorder average, we 
write (JBJ) in the form 

oo 

\nZ L = ^ Wu,i-v(n) W n (u,l-v) + P u ,i-u{n) W n (u,l-u) 

n=0 

+Vv,i-v(n) W n (y,l-v) + zV-uW W n (v, l-u)] (9) 

where f a ,p( n ) is the number of conserving domains of size n bounded by disorder sites at 
densities a and 1 — 13. 

We can average over the configurations of the u aj p(n) by calculating the expectation values 
(v a ^(n)) in the thermodynamic limit (the angled brackets denote a disorder average). This 
is achieved by observing that lini^oo L _1 (z/ Q , i/3 (n)) is just the probability that a site is part 
of an n-site conserving domain bounded by disorder sites with densities a and 1 — (3, hence 

oo 

lim L~ l {\YiZ L ) = (l-p) 2 5> n {q 2 W n (u,l-u)+q(l-q)[W n (u,l-v) 

L— >oo ' * 

n=0 

+W n (v,l-u)] + (l-q) 2 W n (v, 1-v)} (10) 

The form of equation (|10|) is typical of systems which exhibit Griffiths singularities. In 
equilibrium, these singularities are usually inferred from the properties of the Yang-Lee zeros 
of the partition function — we can use the known properties of the Yang-Lee zeros of the 
the analogous quantity, the normalisation Q, to show how Gnffiths singularities anse in 
the disordered nonequilibrium model: For fixed u > 1/2 say, in the complex v -plane and 
for n arbitrarily large, the zeros of (ti|C n |l — v) accumulate arbitrarily close to the point 
v = 1/2 on the real axis. Therefore there exists a singularity in W n (u, 1 — v) arbitrarily 
close to the point v = 1/2 which is exponentially suppressed (by a factor p n ). Thus, such a 
Griffiths-type singularity follows whenever u and v are such that at least one of the W n (a, (3) 
in (jlUJl lies on the phase boundary of the ASEP i.e. whenever u and/or v = 1/2. 

One can go further and consider disorder in the c/s and a/s explicitly. For instance, if 
both Cj and ctj are drawn from binary distributions, then the densities at the disorder sites 
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can assume one of four possible values, each with a different probability in general. However 
Griffiths singularities still arise whenever any one of these values for the density is 1/2, as 
before. 

It is also straightforward to use standard arguments from the study of dilute systems to 
show that the correlation length remains finite at the Griffiths singularity, as is the case in 
equilibrium systems. Note that Griffiths singularities emerge in this non-equilibrium system 
despite the absence of a transition in the pure model. This is in sharp contrast to equilibrium 
systems. 



IV. DYNAMICS: STRETCHED EXPONENTIAL DECAY OF THE DENSITY TO 
ITS STEADY STATE VALUE 

In the pure boundary driven ASEP, whenever the system is in or at the boundary of the 
maximum current phase, the system density decays with time to its stationary value as an 
exponential with a decay constant that depends on system size L as L z , where z = 3/2 is 



the dynamic exponent 
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251 ] . In the low current phase when the boundary injection and 
ejection rates are equal, a shock exists in the steady state, and the dynamic exponent z = 2. 
Otherwise, in the low current phases the relaxation time is finite [24j and does not depend 
on L Hence, in the disordered model, whenever contributions to the normalisation (|1U|) are 
in the maximum current phase, the decay of the system density will be determined by the 
relaxation of these conserving domains. 

For the following analysis, it is sufficient to consider disorder only in the location of the 
disorder sites: we consider homogeneous attachment and detachment rates i.e. Cj = c and 
dj = a. In the case where c = a. (|3*|) gives p = 1/2 at which point conserving domains are 
on the boundary of the maximum current phase; otherwise the conserving domains are in 
the low current phases. Thus only when c = a do we expect the decay constant associated 
with a conserving domain to depend on its sizej^j]. 

Therefore for c = a, in a conserving domain of length n, we assume that the particle 
density p n (t) decays to its steady state value as 

S Pn (t) = p n (t) - 1/2 ~ e- A -\ (11) 

where A n = Aqti^ 3 ^ 2 . In the disordered case, we need to sum over configurations of the 
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disorder sites. This is achieved in the same way as in the previous section so, in the ther- 
modynamic limit, the decay of the system density p(t) becomes 

oo 

5 P (t)~j2p ne ~ Ant > ( 12 ) 

n=0 

where we retain only the n- dependence, as is sufficient to determine the dominant contribu- 
tion to the form of the relaxation. If we convert the sum into an integral and consider late 
times so that the integral can be evaluated at the saddle point, we obtain 

5p(t) ~ exp(-c^), (13) 

where c is a constant and = (1 + z)^ 1 = 2/5. 

Equation (|13p predicts the decay of the density up to some prefactor, power-law in t, 
with an exponent peculiar to the decay of the density. The stretching exponent should 
be universal however, in the sense that other correlation functions e.g. the current, should 
reach their stationary values with the same stretched exponential decay. This result should 
be valid for more general types of disorder whenever one has conserving domains in the 
maximal current phase. 

In figures Q and |2] we show the results of simulations. The simulations were run on 
systems of 10000 sites with periodic boundary conditions and averaged over 1000 histories 
of the dynamics, starting from an empty lattice and with the same realisation of disorder. 
The decay of the averaged system density Sp(t), is shown in Figure [T] for aj = Cj = 10, and 
in Figure El for a,j — c 3 — 1. In both cases the straight line t 2 ^ is given for reference. Figure 
n shows very good agreement with the predicted stretched exponential decay, and even in 
figure |2l where Cj and aj are not large, the agreement is still quite good. Thus it appears 
that our result for the stretching exponent holds for finite rates, although its derivation is 
only exact in the limit of infinite attachment and detachment rates. 

V. CONCLUSION 

In this work we have studied an ASEP with disorder sites where particles are not con- 
served. This may provide a basis for a more realistic model for interacting molecular motors 
such as RNAp. We have used an approximation, exact in the limit of infinite non- conserving 
rates, the underlying assumption behind which is that the system factorises into conserving 
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FIG. 1: Log-log plot of the decay of the density FIG. 2: Log-log plot of the decay of the density 
with time for a 
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domains. According to the ratios of the attachment and detachment rates these domains 
assume different phases of the ASEP with open boundaries. 

Within this approximation we can explicitly identify Griffiths singularities. These arise 
when there are large conserving domains, on the boundary of the maximal current phase. 
An interesting feature is the prediction of a Griffiths singularity despite the absence of a 
transition in the pure system. 

More generally one might ask under what conditions do Griffiths singularities arise in 
nonequilibrium steady states. In equilibrium systems Griffiths singularities are understood 
in terms of Yang-Lee zeros of the partition function. In nonequilibrium systems one does 
not have an energy function, nevertheless one can often identify a quantity that plays the 
role of a partition function, for example the normalisation and recently there has been 
progress in understanding the zeros of such quantities 2^|. 

When there is a spectrum of maximal- current conserving-domain sizes, we have demon- 
strated that correlation functions undergo a stretched exponential decay with a stretching 
exponent predicted to be = 2/5. Moreover simulations suggest this result holds for a wide 
range of attachment and detachment rates. A related stretched exponential decay has al- 
ready been observed for the decay of autocorrelations in a bond diluted symmetric exclusion 
process on a ring (in this case z = 2 so 0=1/3). 

It would be instructive to develop further the approximation that the steady state fac- 
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torises about the disorder sites. As we saw in Section |TT] this approximation is exact in two 
limits, and we gave expressions for the densities at the disorder sites. It would be interesting 
to develop a scheme that interpolates between these two limits. Also of interest would be a 
better understanding of the correlations between the conserving domains which may exist 
away from the two exact limits and their effect on Griffiths singularities. 
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